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Abstract 

The Zipf power law and its connection with the inhomogeneity of the 
system is investigated. We describe the statistical distributions of the 
domain masses in the Potts model near the temperature-induced phase 
transition. We found that the statistical distribution near the critical 
point is described by the power law form with a long tail, while beyond 
the critical point the power law tail is suppressed. 



We use the Potts model [T] for the description of the phase transition. The 
Potts model is the generalization of the Ising model with more than two spin 
components and it has more experimental realizations than the Ising model. 
For a detailed review of Potts model see [2] . 
The Hamiltonian for the q— state Potts model [2] is: 

where Oi G {1, 2, q}, S x>y is the Kronecker delta 



1 if x = y, 
otherwise. 



and 



Jij 



J if i,j are neighbour pairs of spins, 
in opposite case. 



The case q = 2 describes the Ising model. 

The results presented hereafter are obtained by applying the Monte Carlo 
techniques, based on the Metropolis algorithm [3], to the two-dimensional 
q— state Potts model with the periodic boundary conditions and with q = 3 
and 6. 
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A number of exact results for the two -dimensional Potts models are known 
in the infinite volume limit. For example, the phase transition appears at the 
critical temperature T = T c (T c = 2J/k B ln(l + ^Jq) ). It is the second order 
phase transition for q < 4 and the first order one for q > 5, see [I] and [2]. 

The main goal of this paper is the statistical description of the domains in 
Potts model when one approaches the critical point of the phase transition 
induced by the temperature. This problem has been investigated for the Ising 
model in [5]. Our considerations concerning the statistical description of do- 
main masses have universal character and may be used to arbitrary fractal 
system of elements which are described by the random variables x. This vari- 
ables can be listed in the decreasing order. For example: the rank of the city 
connected with its population, the frequency of the occurrence of any word in 
the text, the trading value of largest European's companies or the hyperbolic 
processes in finance [6], [7]. 

We consider the Zipf power law [8] in Bouchaud notation [6]: 

x = k~» (2) 

This law is used often in the description of self-organized critical phenomena. 

In our case x is the number of the Potts spins in the domain called the 
domain mass and k denotes the rank of the domain mass x.The geometri- 
cal clusters are considered [9]. The greatest cluster has the rank 1, smaller 
2 and so on. In the Potts model, clusters are sets of nearest neighbour sites 
occupied by Potts spins with the same orientation. For a given configuration 
of spins clusters are determined uniquely. There is the well establish connec- 
tion between thermal and geometric phase transitions. In the dimension two 
an infinite cluster appears exactly at the critical point. The q— state Potts 
model has the ground state degeneracy q and after a quench from the high- 
temperature phase, small domains start to grow, thus reducing the domain 
boundary curvature. One can notice that there is a difference between the 
domains of the Potts model with q = 2 (the Ising model case) and the Potts 
model with q > 2. For q = 2 the domain boundaries are represented by long 
filaments of one phase within the other one, whereas for q > 2, the domains 
resemble the structure characteristic for polycrystalline grains |10| . [TTj and 
[T2] - i.e. the boundaries are straight and meet at fixed angles. 

We take into account, in our considerations, all domains and we concentrate 
on their size. Bouchaud [B] pointed out the strong correlation between the 
Zipf power law and the inhomogeneity of the system. We are going to test the 
log-log distribution of domain masses versus the rank index k when we come 
near the critical point in the temperature-induced phase transition. 

The slope a of the regression line describing relation between the logarithm 
of the cluster mass and the logarithm of its rank is given by/i ( — = tga), 
and it characterizes the inhomogeneity of the physical structure of the system. 
The inhomogeneity of the system means that its structure has become fractal 
and more hierarchical. We expect that fi value will depend on whether we 
are far or close to the critical point. In particular one should have fi ~ 1 
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in the critical point and fi > 1 beyond the criticality. This follows from the 
following considerations. It is shown in [13] that there is a conjecture between 
Zipf exponent and the Hurst exponent H [14J of the form 

l/fi=\2H-l\ (3) 

It is well known that any complex system exhibits long-range (infinite-long) 
correlations at the critical point. This corresponds to the biggest value of 
Hurst exponent, i.e. H — 1. In such a case one obtains from Eq.(3) fi ~ 1. 
The criticality in the rank statistics was also considered experimentally in 
pJ)],[T6] investigating distributions of island areas of discontinuous metal films 
near the percolation threshold. 

In our consideration we shall concentrate on random variable X which we shall 
call [fi] — variable. We say that the random variable X is [fi] — variable if for 
some xo > and fi > 0, a tail of its distribution function decays as • 

The main property of [p] — variable is that all its moments m q =< x q > 
with q > fi are infinite. The distribution of the cluster size seems to be 
[fi] -distribution and the index fi is a critical exponent. 

The simplest technique which helps to test heavy tail hypothesis and esti- 
mates the tail index fx is based on the following reasoning, (for more details 
see pi]). 

Let X\, X 2 , X n be independent and identically distributed random vari- 
ables with distribution function F. We may re-arrange the Xj in the decreasing 
order. Denoting the largest value, with rank 1, by Xm, the next largest, with 
rank 2, by X( 2 ) etc., we arrive at the sequence of order statistics 

X(i) > X(2) > ... > X( n ). 

Assume that X is a [fi] — variable with tail of its distribution function 
satisfying for some x Q > and fi > 

P{X>x) = (^j\ x>x . (4) 

The assumption (jl]) means that the distribution function of ^- is the Pareto 
distribution with the left endpoint 1 and the parameter fi, and for x > 

P ( uln — > x] = exp(— x) 
\ J 

which means that the distribution function of the random variable a In — is 
exponential with parameter 1. 



3 



Therefore 

Pirn X > x) = exp — — 

V H 

which is a tail of shifted by hi£ and rescaled by /i -1 exponential distribution 
function with parameter 1. 

Drawing the quantile to quantile plot for the exponential distribution func- 
tion we conclude that 



~ ln (jrn7 ' lnXw J ' l - i - n 

should be linear with slope and intercept lnx - In this way we obtain the 
Zipf power law. We may estimate /i in a different way. The most popular 
estimator of fi is the Hill estimator, [18]. For the convenience of the reader 
we recall its definition. The Hill's estimator of - based on m upper-ordered 
statistics > Xp) > ■■■ > X^ m ) is defined as 



H, 



m 



„ = lf>g-^L (5) 



m " X( m+ i) 



Our simulations were performed on two dimensional lattices 600 x 600 with 
periodic boundary conditions. We used the Metropolis algorithm [3]. All 
simulations were initiated with the random initial orientation of spins. We 
thermalized the system to equilibrate the orientations of spins. The number 
of Monte Carlo steps which are necessary to reach the equilibrium state was 
chosen by measuring the average energy of one spin. 

We performed 1 000 000 Monte Carlo steps to equilibrate the system. Once 
the equilibrium of the system is reached we took sample of the cluster con- 
figuration. Then we studied the distribution of the cluster masses. The total 
mass of the cluster was defined as the number of spins in the cluster. 

Our aim is to connect the statistics of domain masses with the way the 
critical point is reached. When we start to advance from the paramagnetic 
phase to the critical point (T — > T c ) (as it is seen in Fig. 1, for q = 3, and 
Fig. 2 for q = 6) the slope of straight lines representing (in log-log scale) the 
Zipf power law increases when fi > 1. 
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Figure 1: The log- log distribution of the domain masses versus the rank order 
index k for the 3-state Potts model for j3 = 0.40 (0), 0.45 (*), 0.49 (□) and critical 
(3 (o), for one configuration. The straight line are with the slopes —1/fi given in 
Table 1. 




12 3 4 

ln(k) 



Figure 2: The log- log distribution of the domain masses Xk versus the rank order 
index k for the 6-state Potts model for (3 = 0.50 (0), 0.55 (*), 0.59 (□) and critical 
(3 (o), for one configuration. The straight line are with the slopes — 1/fi given in 
Table 2. 
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The estimations of indexes -, for the 3-state Potts model with different 
inverse temperature (3 = 1/T, based on 1000 realizations for 600 x 600 lattice 
are presented in Table 1 and Fig. 3. 

Table 1. The estimation of indexes f with standard deviations s. 



p 




s 


0.40 


0.24949 


0.023768 


0.45 


0.30957 


0.030429 


0.49 


0.59078 


0.046805 


Pkr - 0.5025, 


1.09547 


0.047967 




Figure 3: Box plot: The mean value of /i -1 (-) with standard deviation and 95% 
confidence interval, for inverse temperatures (3 = 0.40, 0.45, 0.49 and and close to 
critical 0.50253, obtained by averaging of 1000 Monte Carlo realizations of the 3- 
state Potts model on 600 x 600 lattice. 
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The same estimations of indexes - and its standard deviations as in Table 1 
and Fig 3, but for 6-state Potts model, are presented in Table 2 and Fig.4. 

Table 2. The estimation of index — with standard deviation s. 



p 




s 


0.50 


0.172505 


0.017357 


0.55 


0.193189 


0.020260 


0.59 


0.219017 


0.022300 


f3 kr ~ 0.6191 


0.834464 


0.043940 
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Figure 4: Box plot: The mean value of (-) with standard deviation and 95% 
confidence interval, for inverse temperatures (3 = 0.50, 0.55, 0.59 and close to critical 
0.6191, obtained by averaging of 1000 Monte Carlo realizations of the 6-state Potts 
model on 600 x 600 lattice. 
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Fig. 5 and Fig. 6 represent the dependence between the number of domains 
with mass x normalized by number of configurations and mass(x) for 3-state 
and 6-state Potts models respectively, near the critical point. The distribution 
of domain mass is F[/z]-distribution with /jsil , which means that it satisfies 

power law with a long tail: -T+77, where xq denotes a typical scale and /i ~ 1. 




In(x) 



Figure 5: The dependence between the number of domains with mass x normalized 
by number of configurations (N x /N) and the mass (x) in the log-log scale for the 
3-state Potts model near the critical point T c for 5000 configurations and L = 600. 




Figure 6: The dependence between the number of domains with mass x normalized 
by number of configurations (N x /N) and the mass (x) in the log-log scale for the 
6-state Potts model near the critical point T c for 5000 configuration and L = 600. 



Fig. 7 and Fig. 8 represent the dependence between the number of domains 
with mass x normalized by number of configurations (N x /N) and the mass 
(x) of domain masses for 3-state and 6-state Potts models beyond the critical 
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region when /i > 1, respectively. In this case the distribution of domain masses 
is without heavy power-law tail. This results are in the agreement with the 
standard percolation theory [9]. 
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Figure 7: The dependence between the number of domains with mass x normalized 
by number of configurations (N x /N) and the mass (x) in the log-log scale for the 
3-state Potts model for (3 = 0.4 for 2000 configurations and L = 600. 




Figure 8: The dependence between the number of domains with mass x normalized 
by number of configurations {N x /N) and the mass (x) in the log-log scale for the 
6-state Potts model for (3 = 0.5 for 2000 configuration and L = 600. 



The numerical results presented in Figures 5-8, usually can be deduced from 
the standard percolation theory [9] and are in the agreement with the previous 
results of pj5] and [20]. The origin of the power- law is well understood and 
explained as a consequence of the scale invariance at the criticality regime due 
to the divergence of the correlation length. Beyond the criticality range, where 
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the system is a fractal its correlation length is finite so the power-law tail in 
the distribution is suppressed. One of the possible experimental confirmations 
of the Zipf power law is the distribution of island areas of discontinuous metal 
films, obtained by evaporation of dielectric substances [15], [16]. Another one 
is the distribution of local fields intensities [21] . Unfortunately, the results of 
analogous experiment on ferromagnetic are not known. 

The main result of our paper is obtained by applying the Zipf power law to 
the criticality region and beyond it for the phase transition in the Potts model. 

Further, we have shown, by the simulation results, that the exponent in the 
Zipf power law describes the long tail behaviour of the distribution of domains 
masses at the critical point when /i « 1. 

These results deal with the physical phenomena different from those consid- 
ered e.g., in [13] . [T6] and [21]. However in all cases the formulas describing 
properties of the system are similar and they are based on the Zipf power law. 
The statistical description gives the opportunity to test the phase transition 
not only in the critical region but also beyond the criticality. 
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